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Abstract. In this paper, we introduce local expressions for discrete Mechanics. 
To apply our results simultaneously to several interesting cases, we derive these local 
expressions in the framework of Lie groupoids, following the program proposed by 
Alan Weinstein in [H]. To do this, we will need some results on the geometry of 
Lie groupoids, as, for instance, the construction of symmetric neighborhoods or the 
existence of local bisections. These local descriptions will be particular useful for 
the explicit construction of geometric integrators for mechanical systems (reduced 
or not), in particular, discrete Euler-Lagrange equations, discrete Euler-Poincare 
equations, discrete Lagrange-Poincare equations... The results contained in this 
paper can be considered as a local version of the study that we have started in [T3] , 
on the geometry of discrete Mechanics on Lie groupoids. 



Dedicated to the memory of J. E. Marsden 
1. Introduction 

The use of geometrical methods in the study of dynamical systems (discrete or con- 
tinuous) starts by searching for geometrical structures invariant with respect to the 
given dynamics. It turns out the various such structures emerge naturally for clas- 
sical mechanical systems, as for instance, symplectic or Poisson structures, together 
with various bundle structures. Another geometrical feature that is common to all 
such systems is the presence of symmetries either because there is a redundant or 
extra information in the description of the system or because the system possesses an 
intrinsic invariance. To actually solve them, it is necessary in most occasions to use 
numerical methods. Recently, a new breed of ideas in numerical analysis have come 
that incorporates the geometry of the systems into the analysis and that allows accu- 
rate and robust algorithms with lower spurious effects than the traditional ones (see 
[5] and references therein). Our approach employs the theory of discrete Mechanics 
and variational integrators [12] to derive an integrator for the dynamics preserving 
some of the geometry of the original system. 
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The study of discrete Mechanics on Lie groupoids was proposed by A. Weinstein 
in [121 . This setting is general enough to include discrete counterparts of several 
types of fundamental equations in Mechanics as for instance, standard Euler-Lagrange 
equations for Lagrangians defined on tangent bundles, Euler-Poincare equations for 
Lagrangians defined on Lie algebras, Lagrange-Poincare equations for Lagrangians 
defined on Atiyah bundles, etc. Such discrete counterpart is obtained by discretizing 
the continuous Lagrangian to the corresponding Lie groupoid and then applying a 
discrete variational derivation of the discrete equations of motion. As simple exam- 
ples, for a given different iable manifold Q, the discrete version of the tangent bundle 
TQ is the product manifold Q x Q, equipped with the pair Lie groupoid structure; 
for a given Lie group G, the discrete version of its Lie algebra g is the Lie group G. 

A Lie groupoid G is a natural generalization of the concept of a Lie group, where 
now not all elements are composable. The product gig2 of two elements is only 
defined on the set of composable pairs G2 = {{g, h) E G x G \ (3{g) = a{h)} where 
a : G — > M and (3 : G — > M are the source and target maps over a base manifold 
M. Moreover, in a Lie groupoid we have a set of identities playing a similar role that 
the identity element in group theory. The infinitesimal counterpart of the notion of a 
Lie groupoid is the notion of a Lie algebroid tag '■ — M, in the same way as the 
infinitesimal counterpart of the notion of a Lie group is the notion of a Lie algebra, 
or in other words, the discrete version of a Lie algebroid is a Lie groupoid. 

In [13] we have elucidated the geometrical framework for discrete Mechanics on 
Lie groupoids. In that paper, we found intrinsic expressions for the discrete Euler- 
Lagrange equations, and we have introduced the Poincare-Cartan sections, the dis- 
crete Legendre transformations and the discrete evolution operator in both the La- 
grangian and the Hamiltonian formalism. The notion of regularity has been com- 
pletely characterized and we have proven the symplecticity of the discrete evolution 
operators. The applicability of these developments has been stated in several inter- 
esting examples, in particular for the case of discrete Lagrange-Poincare equations. 
In fact, the general theory of discrete symmetry reduction directly follows from our 
results. 

The main objective of this paper is to obtain local expressions for the different 
objects appearing in discrete Mechanics on Lie groupoids. For this proposal, it is 
necessary to introduce symmetric neighborhoods of a Lie groupoid. A symmetric 
neighborhood is an open neighborhood of one point in the manifold of the identities 
which is "natural" with respect to the structure maps of the Lie groupoid, in the sense 
of Proposition 13.11 Using the coordinates associated to a symmetric neighborhood 
we may write the local expressions of left and right invariant vector fields, ^ and 
^ associated to a section X G T{tag) of the associated Lie algebroid. Now, as we 
have deduced in [9], the discrete Euler-Lagrange equations for a discrete Lagrangian 
Ld-.G^Rare 



Therefore, from this expression, we easily obtain the local expression of the dis- 
crete Euler-Lagrange equations associated to a discrete Lagrangian : G — )■ M, the 
discrete Legendre transformations and we locally characterize the regularity of the 
discrete problem. 
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An interesting point is that when using symmetric neighborhoods we are imphcitly 
assuming that the discrete flow is well defined on this neighborhood. However, this is 
not the more general situation since, in principle, the point Qk and its image Qk+i under 
the discrete fiow may be far enough in such a way both are not included in the same 
symmetric neighborhood. In order to tackle this problem we will use bisections of the 
Lie groupoid which permits to translate neighborhoods of two composable elements 
gk and gk+i to a symmetric neighborhood at the identity point f3{gk) = a{gk+i)- 

The organization of the paper is as follows. In section 2 we recall some constructions 
and results on discrete Mechanics on Lie groupoids which will be used in the next 
sections. In section 3, we will obtain a local expression of the discrete Euler-Lagrange 
equation for a discrete Lagrangian function on a symmetric neighborhood in the Lie 
groupoid where it is defined. In addition we will discuss the existence of local discrete 
Euler-Lagrange evolution operators in such a symmetric neighborhood. Moreover, 
several interesting examples are considered. The existence of general local discrete 
Euler-Lagrange evolution operators is studied in section 4. For this purpose, we 
will use bisections on the Lie groupoid. The paper ends with our conclusions and a 
description of future research directions. 



2. Groupoids and discrete Mechanics 

2.1. Lie groupoids. In this Section, we will recall the definition of a Lie groupoid 
and some generalities about them are explained (for more details, see [3l [T2]). 
A groupoid over a set M is a set G together with the following structural maps: 

• A pair of maps a : G ^ M, the source, and /3 : G — )■ M, the target. Thus, 
an element g & G is thought as an arrow from x = a{g) to y = /3{g) in M 



x = a(g) y = l3{g) 

The maps a and /3 define the set of composable pairs 

G2 = {{gM^GxG/l3{g) = a{h)}. 

A multiplication m : G2 ^ G, to be denoted simply by m{g, h) = gh, such 
that 

— a{gh) = a{g) and f3{gh) = (3{h). 

- gihk) = {gh)k. 

If g is an arrow from x = a{g) to y = (5{g) = a{h) and h is an arrow from y 
to z = (3{h) then gh is the composite arrow from x to z 




' 0(9) = a{gh) 



y = P(g) = a(h) 



z = P{h) = p{gh) 



An identity map e : M G, a. section of a and /3, such that 

- £ioi{g))g = g and ge{(3{g)) = g. 

An inversion map i : G ^ G, to he denoted simply by i{g) 
that 

- 9~^9 = ^Wig)) and gg'^ = e{a{g)). 



such 
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A groupoid G over a set M will be denoted simply by the symbol G ^ M . 

The groupoid G ^ M is said to be a Lie groupoid if G and M are manifolds 
and all the structural maps are different iable with a and (3 differentiable submersions. 
If G =^ M is a Lie groupoid then m is a submersion, e is an immersion and z is a 
diffeomorphism. Moreover, if a; G M, a~^{x) (resp., (3~^{x)) will be said the a-fiber 
(resp., the (5-fiber) of x. 

On the other hand, g E G then the left-translation by g E G and the right- 
translation by g are the diffeomorphisms 

Ig : a-\l3{g)) — > a~\a{g)) ; h — > lg{h) = gh, 
rg : l3-\a{g)) ^ l3-\/3{g)) ; h^rg{h) = hg. 

Note that = Ig-i and = Vg-i. 

A vector field X on G is said to be left- invariant (resp., right- invariant) if 
it is tangent to the fibers of a (resp., /3) and X{gh) = (Thlg){Xh) (resp., X{gh) = 
{Tgrh){X{g))),foT{g,h)eG,. 

Now, we will recall the definition of the Lie algebroid associated with G. 

We consider the vector bundle r : AG M, whose fiber at a point x G M is 
AxG = Ve(x)Ci = Ker{T^(^x)Oi). It is easy to prove that there exists a bijection between 
the space r(r) and the set of left-invariant (resp., right-invariant) vector fields on G. If 
X is a section of r : AG — )■ M, the corresponding left-invariant (resp., right-invariant) 

vector field on G will be denoted ^ (resp., ^), where 

(1) ^(^?) = {TeiPig))lg){X{m)), 

(2) '^{g) = -(r,(„(,))r,)((T,(,(,))z)(X(a(^?)))), 

for g E G. Using the above facts, we may introduce a Lie algebroid structure ([■,■], p) 
on AG, which is defined by 

(3) txyj = H P{X){x) = (T,(,.)/3)(X(x)), 
for X,Y e r(r) and x E M. Note that 

(4) lX^^ = -[?t,^], [^,^] = 0, 

(5) TioJ^ = oz, Tio'X = o z, 
(for more details, see [H fT2]). 

2.2. Discrete Euler-Lagrange equations. Let G be a Lie groupoid with structural 
maps 

a,/3:G^M, e:M^G, i:G^G, m : G2 ^ G. 
Denote by r : AG — j- M the Lie algebroid of G. 
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A discrete Lagrangian is a function L^: G — M. Fixed g E G, we define the set 
of admissible sequences with values in G: 

= {{g,, . . . ,^?^) G / {gk,gk+i) eG2foTk = l,...,N-l 

and gi...gn = g}. 

An admissible sequence {gi, . . . , g^) G is a solution of the discrete Euler- 
Lagrange equations if 



N~l 

, , ^ ,, ,^ forXi,...,X^_iGr(r). 

k=l 



For = 2 we obtain that ((7, h) G G2 is a solution if 

1l{g){Ld)-l{h){Ld)=^ 

for every section X of AG. 

2.3. Discrete Poincare-Cartan sections. Given a Lagrangian function L^: G 
M, we will study the geometrical properties of the discrete Euler-Lagrange equations. 
Consider the vector bundle 

tt" : P^G = Vl3®Va ^G 

where Va (respectively, V/S) is the vertical bundle of the source map a : G — >■ M 
(respectively, the target map (3 : G ^ M). Then, one may introduce a Lie algebroid 
structure on vr^ : P^G = VP®Va^G. The anchor map p^"^ : P^G = Vf3®Va^ 
TG is given by 

p^^«(X„ = X, + Yg, for (X„ Yg) G T/^/? © K,a 

and the Lie bracket |-, ■]^^'^ on the space r(7r^) is characterized by the following 
relation 

(6) [(^, H (x\ hr^ = i-ix,x'i fi^). 



for X,Y,X',Y' G T{t) (see [I3]). 

Now, define the Poincare-Cartan 1-sections 0^^,6^^ G r((7r^)*) as follows 

(7) Ql^{g){Xg,Yg) = -X,(L,), e+^(^?)(X„ F,) = F,(L,), 

for each 5^ G G and (X^, F^) G V^/? © VgU. 

If (i is the differential of the Lie algebroid vr^ : P^G = V [i © — )■ G we have 
that dLd = Qt — ©7 and so, using (f = 0, it follows that dQt = dQ7 . This 

" ' ^ ' ^d id 

means that there exists a unique 2-section fi^^ = —dQ\^ = —dQ^^y that will be 
called the Poincare- Cartan 2-section. This 2-section will be important for studying 
symplecticity of the discrete Euler-Lagrange equations. 

Let X be a section of the Lie algebroid r : AG — )■ M . Then, one may consider the 
sections X^^'^^ and X^^'^^ of the vector bundle vr^ : P^G = V (3 ®Va ^ G given by 

X(i'°)(^7) = (^(^7),0,), X(°'^)(^7) = (0„^(^?)), for ^?gG. 
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Moreover, if (7 G G, {X^} (respectively, {V^}) is a local basis of r(r) in an open 
subset U (respectively, V) of M such that a{g) G U (respectively, (3{g) G V) then 
{x!^^'^\Yi^'^^} is a local basis of r(7r") in a-\U) n ^'^{V) and 

(8) fiL,(x^^ = i^..(xf ^\ = 0, 

and 

(9) ^lAX^'''\y('''^) = %{X!,{L,)) = X'.fyl,)). 
(for more details, see [I3]). 

2.4. Discrete Lagrangian evolution operator. We say that a differentiable map- 
ping \E' : G — y G is a discrete flow or a discrete Lagrangian evolution oper- 
ator for Ld if it verifies the following properties: 

- graph(^^) C G2, that is, (^f, '^{g)) G G2, V^f G G. 

- [g, "^{g)) is a solution of the discrete Euler-Lagrange equations, for all g E G, 
that is, 

(10) ^(^?)(L,)-^(v[/(^))(L,) = 
for every section X of AG and every g E G. 

2.5. Discrete Legendre transformations. Given a discrete Lagrangian L^: G ^ 
M we define two discrete Legendre transformations ¥~Ld : G — yl*G and 
¥+Ld : G A*G as follows (see [13j) 

(11) (¥'Ld)ih){Ve{a{h))) = -Ve{a{h)){Ld O O i) , for Ve{a{h)) ^ ^a(/i)G, 

(12) (¥+Ld){g){v,^,3(g))) = v^^p(g^){Ld o Ig), for 17^(^3(3)) G A;3(g)G. 

Remark 2.1. Note that (F+Ld)(^) G A^(g)G and (F-Lrf)(/i) G A^^^^G. Furthermore, 
if {X^} (respectively, {1^^}) is a local basis of r(r) in an open subset U such that 
a{h) G U (respectively, I3{g) G V) and {^'^} (respectively, {i^^}) is the dual basis of 
r(r*), it follows that 

¥-Ld{h) = jl^{h){Ld)X^{a{h)), ¥-^Ld{g) = ^ ^{g){Ld)Y>^{f3{g)). 

o 

2.6. Discrete regular Lagrangians. A Lagrangian : G — )■ M on a Lie groupoid 
G is said to be regular if the Poincare-Cartan 2-section is symplectic on the Lie 
algebroid vr"^ : P'^G = V/3 (Bg Vc( G, that is, is nondegenerate (see [TB]). 

Using ([H]), we deduce that the Lagrangian Ld is regular if and only if for every 
g E G and every local basis {-^^7} (respectively, {V^}) of r(r) on an open subset U 
(respectively, V) of M such that a{g) G U (respectively, (3{g) G V) we have that the 

matrix X^(l^(Ld)) is regular on a~^{U) fl /3~-^(V^). 

In [13], we have proved that the following conditions are equivalent: 

• Lrf : G — > M is a regular discrete Lagrangian function. 

• The Legendre transformation ¥~Ld is a local diffeomorphism. 

• The Legendre transformation F+L^ is a local diffeomorphism. 
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Moreover, if : G — j- M is regular and {qq, Hq) G G2 is a solution of the discrete 
Euler-Lagrange equations for then there exist two open subsets Uq and Vq of G, 
with G Uq and Hq G Vq, and there exists a (local) discrete Lagrangian evolution 
operator "^^^ : Uq ^ Vq such that: 

• *Ld(^o) = ho, 

• is a diffeomorphism and 

• is unique, that is, if Uq is an open subset of G, with G Uq and : 
t/g — )■ G is a (local) discrete Lagrangian evolution operator then "^'^ IUqCiUq = 

3. Discrete Euler-Lagrange equations: symmetric neighborhoods 
3.L Symmetric neighborhoods. First, we prove the following result 

Proposition 3.1. Let lA he an open subset of G and xq E M he a point such that 
e{xo) G U. There exists an open suhset W CU of G with e{xo) G W and such that 

(1) e{a{W)) C W and e(/3(W)) C W, 

(2) i(W) = W, and 

(3) m((W X W) nCa) CU. 

The open suhset W is said to he a symmetric neighborhood associated to U and 

Xq. 

Proof. The multiplication map m: G2 G is continuous, so that we may choose 
an open subset Wi of G such that e^xq) G Wi and m((Wi x Wi) fl G2) C U. Since 
the identity map e: M — )■ G is also continuous, we deduce that there exists an open 
subset V of M such that xq G V and e{V') C Wi. Thus, if we consider the open 
W2 of G given by W2 = Wi n a~^(V') n /3~HV'), then it is clear that e{xo) G W2 
and moreover it is easy to prove that e(a(W2)) C W2 and e{f3{'W2)) C W2, and also 
m((W2 X W2) n G2) C U. Finally if we take W = W2 n ^(^2) it follows that >V 
satisfies the three above mentioned conditions. □ 

3.2. Local coordinate expressions of structural maps. On a symmetric neigh- 
borhood of a point it is easy to get local coordinate expressions for the structure maps 
of the Lie groupoid G. We consider a point Xq G M and a local coordinate system 
(x, u), defined in a neighborhood W C G of e{xo), adapted to the fibration a : G — )■ M, 
i.e. if the coordinates oi g ElA are (x', u'^) then the coordinates of a{g) G M are (x*). 
We can moreover assume that the identities correspond to elements with coordinates 
(x,0). The target map /3 defines a local function b as follows: if the coordinates of 
g are (x,u), then the coordinates of I3{g) are b(x, u). Note that b(x,0) = x. Two 
elements g and h with coordinates (x, u) and (y, f ) are composable if and only if 
y = b(x,M). Hence local coordinates for G2 are given by {x,u,v). 

To obtain a local description for the product, we consider a symmetric neighborhood 
W associated to xq and U. If two elements g,h eW with coordinates (x, u) and {y, v) 
respectively, are composable then y = b(x,M), and the product gh has coordinates 
(x, p(x, M,f)) for some smooth function p. We will write 



(13) 



(x,u) ■ {y,v) = (x, p{x,u,v)). 
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The relation f3{gh) = /3{h), for {g, h) G 6*2, imposes the restriction b(?/, v) = b{x, p(x, u, 
i.e. 

(14) b{b{x, u), v) = b(x, p(x, u, v)). 

The property ge{/3{g)) = g, for g G W, is locally equivalent to the equation p(x, u, 0) = 
u] while the property e{a{g))g = g is locally equivalent to the equation p{x, 0, v) = v. 
Therefore 

(15) p(x, M,0) = M, p{x,0,v) = v. 
Associativity {gh)k = g{hk) imposes the further relation 

(16) p{x,p{x,u,v),w) = p{x,u,p{y,v,w)) with ?/ = b(x, m). 

In what follows we will use the following functions defined in terms of b(x, u) and 
p{x,u,v), 

du'y' 



(17) Llix,u) = ^ix,u,0) 



dpf 

We will also take into account that 

^ ' dp^ 92 p^ 

which follow from (|T5|) . The only relevant second order derivatives are given by 

(19) ClXx) = ^ (x, 0, 0) - ^ ; (x, 0, 0). 

From the definition of L^! and Kj, it follows that 

.f2,.,o)-|??(..o,. 

On the other hand, if i : G — )■ G is the inversion we have that 

i{x,u) = {b{x,u), l{x,u)) 
and the condition i{e{x)) = e{x), for all x G M, implies that 

l{x,0) = 0. 

Moreover, using that p(a;, u, l{x, u)) = 0, we deduce that 

dp^ di" dp^ 

^{x, u, l{x, u)) + ^{x, u)^{x, u, l{x, u)) = 0. 

Thus, from f lTSjl . we obtain that 

(21) ^^{x,0) = -6;. 
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3.3. Invariant vector fields. The local expression for left- and right-translations are 
easy to obtain. For go & W C G the left translation Ig^ is the map Ig^ : a~^{/3{go)) — )■ 
a~^{a{go)), given by Ig^^g = g^g. If g^ has coordinates (xo,Mo), then the elements on 
the a-fiber a~^{l3{gQ)) have coordinates of the form (b(xo, Mq), "y), and the coordinates 
of Igf^g are (xq, p(a;o, uq, v)). We will write 



(22) 



kxo,uo){Kxo^Uo),v) = (Xo, p(Xo,Mo,w))- 



Similarly, for /iq G W C G the right translation map rh^: f3~^{a{hQ)) — )■ (3~^ {(3 (hQ)) , 
is defined by Vh^g = gh^. If Hq has coordinates {xo,uq), then the elements on the /3- 
fiber /3~^(a(/io)) have coordinates of the form [x, u) with the restriction b(x, u) = xq, 
and the coordinates of Vh^g are (x, p(x,m,mo)). We will write 

(23) r(^,^^uo){.x, u) = (x, p(x, u, Uq)). 

A left-invariant vector field is of the form ^{g) = T^(i3[g))lg{v) for v G ker T£(^(g))Q;. 
To obtain a local basis of left-invariant vector fields we can take the local coordinate 
basis = ■^\s{i3{g)) of ker T£(^(g))a. Thus, ioi g E G with coordinates (x, u), we have 



(24) t-,{g) 



d 



9p^ 9 

(X,M,0)- 



{x,u) 



Similarly, a right-invariant vector field can be written in the form ^((7) = T^[oi{g))f^ g{v) 
for V G ker T£(Q,(g))/9. To obtain a local basis of right-invariant vector fields we first 
have to look for a basis of the vector space ker Te(Q,(g))/3. From the definition of the 
functions p!j,, it follows easily that the vectors = Pj^ — -^ are in ker Te(Q,(g))/3, and 
moreover they are related to the vectors by the inversion map, that is Ti{e^) = f^. 
It follows that a basis of right invariant vector fields is given by 



^{9) = Teia{g))rg { -p\^ 



( ^ ^ 




d 




V ^''dx' 


H 

e(c^(9)) 


du-y 


£("(9))) 



(25) 



_d_ 
(9x* 
_d_ 
(9x* 



-pl{x)——{x,0,u) + TT— (a:, 0,m' 



dx^ 
d 



du-y 



d 



g dul" 

where as before (x, u) are the coordinates for g E G. Note that, from f|T5|) . we deduce 
0. 



dp^" 

that — — (x, 0, u) 



dx^ 



3.4. The Lie algebroid of G. The Lie algebroid of G is defined on the vector 
bundle t: E ^ M with fiber at the point x G M given by = kerT£(j.)a. A local 
basis of sections of E is given by the coordinate vector fields e^(x) = ^ 



anchor is the map p: E 
local coordinates, if a = 



TM defined by p(a) 



then p{a) = pl^{x)y'y-£. 



^\e{x)- The 
^£(x)/3(fl), where x = r(a). In 

The bracket is defined in 



terms of the bracket of left- invariant vector fields. A simple calculation shows that 
[^, t^] = G^^^i,, with C^^ given by ([H]), from where we get [e-^, e^] 



10 



J. C. MARRERO, D. MARTIN DE DIEGO, AND E. MARTINEZ 



3.5. Discrete Euler-Lagrange equations. Consider now a discrete Lagrangian 
function L^. A composable pair {g, h) G G2 satisfies tfie Euler-Lagrange equations 
for Ld if 

(26) ^{g){Ld) = '^{h){Ld) for every section X of E. 

If botli g and h are on the same symmetric neighborhood W, with coordinates (x, u) 
for g and {y, v) for /i, we can apply the local results above and we readily get the 
coordinate expression of the Euler-Lagrange equations 

f)T r)T f)T 

(27) L^,ix,u)^ix,u)+pl^{y)^{y,v) - R',{y.v)^^{y.v) = 

(28) y = b{x,u), 

where the second equation takes into account that P{g) = a{h). 

Assume that we have a solution {go, ho) e G2 of the Euler-Lagrange equations. To 
analyze the existence of solution of the Euler-Lagrange equations for elements {g.,h) G 
G2 near {go, ho), we can apply the implicit function theorem. In the application of 
such a theorem, the relevant matrix [(FLrf)^]^^^=i m is the following 

(29) iWLd)lM = p;(-)^(-,-) - 

Proposition 3.2. Let {yo,VQ) be the coordinates of the point Hq. The following state- 
ments are equivalent. 

• The matrix {¥Ld)'p{yo,VQ) is regular. 

• The Poincare-Cartan 2-section Ql^ is non-degenerate at the point Hq. 

• The map ¥~Ld is a local diffeomorphism at Hq. 

Any of them implies the following: There exist open neighborhoods Xq C W and 
yo C W of go and ho such that if g ^ Xq then there is a unique "^{g) = h & yo 
satisfying that the pair {g, h) is a solution of the Euler-Lagrange equations for Ld- In 
fact, the map : Xo ^ yo is a local discrete Euler-Lagrange evolution operator. 

Proof. We first notice that the local expression of the map F^^ is 

^{x,u) + R^^{x,u)—. 

The differential of this local function at the point h = {yo, Vq) is of the form 



-¥Ld{yo,vo) 

from where the equivalence of the first and the third assertions immediately follows. 

For the second we just take a local basis of sections {e^} of E, defined in a neigh- 
borhood of Xq and associated to the a- vertical vector fields and we compute the 

value of on the associated basis {eiy'^\e^'^^}. From and it follows that 
^Ld(^o) is regular if and only if the matrix 

n,,iyo,Vo) = fiz.,(e(^'°\e!f'i))(/io) = ^(4(i:,))(/io) 



(F^J(x,m) = ( x,-pl^{x)——{x,u) + R':;{x,u)^{x,u] 
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is regular. From the expressions (1241) and fl25l) this matrix is 



= -Lj(yo,t^o)(FL,)^(2/o,^^o). 

It follows from ^IM that the matrix L^^{yo,vo) is regular. Thus, we get that the 
regularity of the matrix (¥Ld)^{yQ,vo) is equivalent to the regularity of Qi^{ho). 

Finally, let be the left-hand side of the discrete Euler-Lagrange equations ({2] 
once the equation (l28i) has been used 

dL dL 
(30) A^(x, u, v) = Ll{x, u)^{x, u) + p;(b(x, u))^{b{x, u),v) + 



Rl{h{x,u),v)^{h{x,u),v). 



did 

We want to study the existence of solution of the discrete Euler-Lagrange equations 
\^{x,u,v) = in a neighborhood of the point {xq.uq^vq), where (x, m) are the data 
and V is the unknown. Applying the implicit function theorem we have to study the 
regularity of the matrix ^^(xo, Uq, Vq). A straightforward calculation shows that this 
matrix is just 

d\ 

^(xo, uo, ^'o) = {¥Ld)liyo, vq) with yo = b(xo, uq), 

from where our last assertion readily follows. □ 

Remark 3.3. Suppose that G e{M). Then, ho has local coordinates {yo,0). 
Moreover, from flTTI) and ffTSl) . it follows that the matrix L^(|/o, 0) is the identity 
matrix. Therefore, we deduce that the matrix Q^^{yo, 0) is, up to the sign, the matrix 
(¥Ld)ii^{y 0,0) (see the proof of Theorem 14.91 in Section l42l) . o 

In general the Euler-Lagrange equations are understood as the equations deter- 
mining h = (yiv) from the already known data g = {x,u), as we did in the proof of 
the above theorem. However, we can also try to solve these equations backwards to 
obtain {x,u) from {y,v). Instead of applying again the implicit function theorem to 
the system of equations (1271) and (125]) we can rewrite the equations in a different co- 
ordinate system, adapted to the fibration /3, as follows. On the open subset U = i(U) 
we consider the local coordinates {x,u) defined as {x,u) = {x,u) o z, where i is the 
inversion map in the groupoid. Since P o i = a, we have that these new coordinates 
are adapted to the submersion /3 and we can proceed as in the previous case. In these 
coordinates, a basis of the left-invariant and right-invariant vector fields is 



(31) t;(,,^^-p-(,)^^^^ + i^(,,^g^ 

d 



d 



{x,u) 



(32) = L^(,,^^^ 

3.6. Examples. We present here some illustrative examples. 
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3.6.1. Pair or Banal groupoid. We consider as a first example the pair (banal) groupoid 
G = M X M, where the structural maps are 

a{x,y) = x, I3{x,y) = y, e{x) = {x,x), i{x,y) = {y, x), 
miix^y), {y,z)) = {x,z). 

The Lie algebroid of G is isomorphic to the standard Lie algebroid tm '■ TM — )■ M, 
therefore the pair groupoid is considered as the discrete phase space for discretization 
of Lagrangian functions L : TM — )■ M. 

Let Xo G M and a local coordinate system (x*) defined on a neighborhood V' of 
xq. Then 14 = V x V' is obviously a symmetric neighborhood. For a fixed h > 0, an 
associated coordinate system adapted to the a-projection is 

x^{xi,X2) = x'^{xi) and u\xi,X2) = • 

Obviously the identities correspond to the elements (x*,0) in this coordinate system 
and a{x,u) = x and P{x,u) = x + hu. On the other hand, the composition of two 
elements (x, m) and {x + hu,v) is {x,u + v). Therefore p(x, = u + v. The inversion 
map is now i{x, u) = {x + hu, —u). 

If we take the natural local coordinate basis {g^|e(a;i)} of kerTe(2;i)a then 

V _ 

d ^ d d 
du^ dx'- (9m* 
Therefore, the discrete Euler-Lagrange equations are 

^(x,n) + /.^(y,.)-^(,,.) = 
y = X + hu. 

The discrete Lagrangian is regular if the matrix 



dx^du^ du^du^ 
is non singular. 

Example 3.4. As a concrete simple example, consider the continuous Lagrangian 

L{x,x) = —x'^Mx — V{x) 

(with M a constant symmetric invertible matrix). A typical discretization for the 
lagrangian is, for instance, 

Ld{x, u) = ^u^Mu — hV{x + "^^)- 

Then, the discrete Euler-Lagrange equations are: 

,,v-u 1 fdV , h ^ dV , h , 

y = X + hu, 
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which leads us the classical implicit midpoint rule. 
If on the other hand, we take the discretization 

Ld{x,u) = ^{L{x,u) + L{x + hu,u)) 

= ^u^Mu - ^ {y{x) + V{x + hu)) , 

then the corresponding discrete Euler-Lagrange equations are: 

^^v — u dV , . 

y = X + hu, 

which is a representation of the Stormer-Verlet method. 

3.6.2. Lie groups. Another interesting example corresponds to the case of Lie groups. 
In this case, we consider a Lie group G as a groupoid over one point M = {e}, the 
identity element of G. The structural maps are 

= 1^(9) = ^, £^(e) = e, ^Ifi-) = 'm{g,h) = gh, ioi g,heG. 

The Lie algebroid associated with G is just the Lie algebra q = T^G of G. 

Near of the identity, it is interesting to regard the elements G G as small dis- 
placements on the Lie group. Thus, it is possible to express each term through a Lie 
algebra element that can be regarded as the averaged velocity of this displacement. 
This is typically accomplished using a retraction map r : g — )■ G which is an analytic 
local diffeomorphism around the identity such that T(^)r(— ^) = e, where ^ G (see 
[2]). Thereby r provides a local chart on the Lie group. 

Given a retraction map r its right trivialized tangent dr^ ^ — and is inverse 
dr^"^ : — ^ are defined for all 77 G g as follows 

The retraction map allows us to transport locally the Lie group structure on an 
open symmetric neighborhood of G to a local r-dependent Lie group structure on the 
Lie algebra g defined on a local neighborhood ^ of G g. We will write T{h^) = g 
for an enough small time step h > such that h^ & V q. 

Now, on the Lie algebra it is easy to consider local coordinates because it is a vector 
space. In consequence, fixing a basis {e^} of g, we induce coordinates {u'^) on 0. 

In these coordinates, a basis of left- and right-invariant vector fields is 

^iv) = Tr{hr,)r~\T{hr])e^) = dTf;^\Adr(hv)e^) 
e^iv) = T^ihrj)T^\e^T{hr])) = dr^^\e^), 

where ?7 G g. Given a lagrangian / : g — )■ M, we deduce that the discrete Euler-Poincare 
equations are: 

Kivk){i) - 4(%+i)(0 = 0, 

or, alternatively, as discrete Lie-Poisson equations: 
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dl 

Example 3.5. As an example of retraction map we will consider the Cayley transform 
which is one of the most computationally efficient parametrizations of a Lie group. 
The Cayley map cay : — G is defined by cay(,^) = (e — |)~^(c + |) and is valid for 
a general class of quadratic groups as for instance SO{n), SE{n) or Sp{n) (see [5]). 
Its right trivialized derivative and inverse are given by 

dcay^r/ = (e - (e + ^)"^ 

dcay^S = (e- -)r/(e + -). 

In this case it is possible to write more explicitly the Lie-Poincare equations. In 
fact, we have that 



5tM = (<-'^)eA'+'^: 



with 7] G 50 (3). Therefore, the discrete Euler-Poincare equations are: 



h fi^ dl 

= (^7 - 2 [^^' ~ -jVke^iVk, Q^iVk)) 

h h"^ dl 

-(^7 + 2^(^i^Vk+i] - — r/fc+ie^r/fc+i,— (r/fe+i)). 

As a concrete example, we consider the case of the group 5*0(3), using the Cayley 
transformation and fixing the standard basis of its Lie algebra so (3): 



/OO 0\ /00l\ /0-10\ 

ei = -1 , 62 = , 63 = 1 
\01 0/ \-100/ \0 00/ 
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Thus, the elements r]k, Vk+i in will have coordinates (x^, i/k, Zk) and (x^+i, Vk+i) ^fc+i), 
respectively, in this fixed basis. This allows us to write the previous discrete Euler- 
Poincare equations as follows: 







1 + ^1^ 



+ 



+ 



dL 

dx 

hyk _ h'^XkzA dL 
2 4 J dz 

hzk+i h'^Xk+iVk 



4 



hzk 
2 



1 + 



+1 



+ 



2 

ffxkVk 
4 

hxk 



h'^VkZk 



hyk 
2 



Uk+l 
4 

h^XkZk 



4 

dL 

dx 

dL 

dz 

dL 

dy 
dL 
dx 



dL 
dy 



k+l 



\ dL 




J dy 


k 


^) 


dL 




dx 


f hyk+i 


\ 2 



k+l 
2 



1 



hWk 



dL 
dy 



h^Xk+iZk+i 
4 



dL 

dz 



k+l 



k+l 



hzk+1 h?Xk+iyk+i 
2 4 

hxk+i 



h'^yk+iZk+i 



dL 

dx 

dL 

dz 



k+l 



k+l 



hxi 



h'^ykZk 



h'^Zk 



dL 

dz 



where L 



hxk+i K^yk+iZk+i 
2 4 

is defined by L{x,y,z 



2 

hyk+i 

2 

dL 
dy 



dL 

dy 



h'^Xk+iZk+i 



k 

dL 

dx 



k+l 



1 + 



k+l 

dL 

dz 



k+l 



l{xei + ye2 + ze-^] 



3.7. Transformation or action Lie groupoid. The theory of Lie groupoids covers 
other interesting examples that are useful for the construction of variational integra- 
tors for different mechanical systems. This is the case of transformation or action Lie 
groupoids. Let G be a Lie group and ■ : M x G — )■ M, (x, ^) G M x G H- xg, a right 
action of G on M. Consider the Lie groupoid G = M x G over M with structural 
maps given by 



(33) 



a{x,g)^= x,^ (3{x,g) 
m{{x,g), {xg,g')) = ( 



= xg, 

x,gg'), 



e{x) = {x, e), 
and i{x,g) = {xg,g-'^). 



where e is the identity on G. The Lie groupoid G is called the action or transformation 
Lie groupoid. Its associated Lie algebroid is the action algebroid pri : M x g — > M 
where g is the Lie algebra of Lie group G (see [13] for details about the Lie algebroid 
structure). We have that r(pri) = {f/ : M — )■ g | where f] is smooth mapping}. In 
particular, if 77 G g then rj defines a constant section G,, : M — g of pri : M x g — )■ M: 
Gri{x) = {x,ri) for all x G M. It is possible to check that the corresponding left- and 
right-invariant vector fields on G are (see [T3]): 

(34) £^,(x,^) = (0.,V(^)), ^ 



X, 



-Vm[x 
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for {x,g) E G = M X G and where tjm is the infinitesimal generator of the right action 
■ : M X G — )■ M associated with 77. 

Let Ld : G = M X G ^ M. he a. discrete Lagrangian. Then, a composable pair 
{{xk,gk), {xkgk, gk+i)) G 6*2 is a solution of the discrete Euler-Lagrange equations for 



Given a retraction map r : g — t- G, we can transport the Lie group structure on 
an open symmetric neighborhood of e in G to an open neighborhood of in g as in 
Subsection 13.6.21 Moreover, if xo G M, we may assume, without loss of generality, 
that this neighborhood acts on an open neighborhood U of the point Xq G M. Thus 
if ^ G the fiow of the fundamental vector field C,m associated with this local 
action is given by 



for t G M and x E U. Therefore, as in the previous section, fixed a basis {e^} of q, 
we induce a basis of left- and right-invariant vector fields 



where f] E g and x E U. Given a lagrangian Z : M x g — t- M, we deduce that the 
discrete Euler-Lagrange equations are: 

= dr^";^t(Adr(fc*)fe)eT,)(/^J - dT^l^^{e^){l^^^J + h{e^)M{xk+i){lfj^^,), 

Xk+l = XkT{hflk) . 

where for every f/ G g (resp., x G M) we denote by If, (resp., Ix) the real function on 
M (resp., on g) given by lf,{y) = l{y,fj) (resp., lx{fi') = l{x,fj')). 

Example 3.6. As a typical example of a discrete system defined on a transformation 
Lie groupoid consider a discretization of the heavy top [TOl [T6l [13]. This system is 
modelled on the transformation Lie algebroid r : S*^ x so (3) — ?■ S"^ with Lagrangian 



where G ~ so (3) is the angular velocity, F is the direction opposite to the 
gravity and e is a unit vector in the direction from the fixed point to the center of 
mass, all them expressed in a frame fixed to the body. The constants m, g and d are 
respectively the mass of the body, the strength of the gravitational acceleration and 
the distance from the fixed point to the center of mass. The matrix I is the inertia 
tensor of the body. 

We will also use the Cayley transformation on 5*0(3) to describe the discrete Euler- 
Lagrange equations for the heavy top. We have that 



Ld if 




<^,Jt,x)=xexpa{t{hO) 
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with 



/ -ils 1^2 \ 

r]= { -1^1 , 

\ -^2 y 



{e^} the standard basis on S0{3) and T G S'^. 

Therefore, the discrete Euler- Lagrange equations are: 



" ' (^fc+l)) 



(e + 2 



9r 



(r 



Tfc+i = Tfe cay (/ir/fc), 



or, in other terms 



. ,,(n,uii + ^Vm)/^^<"'W"=>'.'""''^ 



2jA:+l 



^(^3)fe+l /i^(^l)fc+l(^2)fc+ 



3jfc+l 



2 4 

/l(l]2)fc+l , /l2(fil)fc+l(fi3) 



2 

+hmgd{Zk+ie2 - Ifc+iCg) 



r ^ ^ ^^(^2)fc(^^3)A , ,o ^ fh^k+i , /i'(l^i)fc+i(fi2)fc+i 



2 47' V 2 4 

-'2l"2jfe+l J- H ^ + -'3l"3jfc+l 



4 

-hmgd{Xk+ie-i - Zt+iei) 
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) 



-hmgd{Xk+ie2 - Ifc+iei) 



= (Xfe+i,n+i,Zfc+i)-(X,,n,Z,)(e-^)"i(e + ^) 



where Tk = (X^, Yk, Zk) G with X| + + Zl = l,e = (ei, ea, es) and 



/ -(l]3)fc {^2)k \ 

Vk = i^3)k 

\ -{n2)k {^i)k / 



Remark 3.7. Our approach also admits other interesting examples. For instance, 
assume that we have a discrete system modeled hj a. : M x M x G ^ K. which is an 
approximation of a continuous Lagrangian L : TM x g — R. This lagrangian typically 
appears as reduction of a G- invariant Lagrangian function L : T{M x G) — R (see 
[13], for the general case). Of course we can combine the techniques exposed in 
Subsections 13.6.11 and I3.6l2l to obtain a local description of the corresponding discrete 
Euler-Lagrange equations in terms of the continuous Lagrangian L. o 

4. Bisections and discrete Euler-Lagrange evolution operators 

One of the main limitations of the techniques exposed in Section |3] is that we need 
to work in a neighborhood of the identities of the Lie groupoid G. Using an enough 
small time stepping we can guarantee that the evolution of the evolution operator for 
a discrete Lagrangian takes values on the chosen symmetric neighborhood, even it is 
possible to adapt the time stepping to make it happen. Another possibility is to use 
the notion of bisections on Lie groupoids. As we will see it will allow us to consider 
points far from the identities completing our local description of discrete Mechanics. 

We consider now the general case of a solution [go, ho) G G2 of the Euler-Lagrange 
equations where the points Qq and Hq E G are not necessarily close enough to be con- 
tained in a common symmetric neighborhood. If we want to obtain a local expression 
of the discrete Euler-Lagrange operator which connects the above points, we must 
choose suitable neighborhoods of go and Hq. For this purpose, we will consider a sym- 
metric neighborhood W, a local bisection through the point go and a local bisection 
through ho- By left-translation and right-translation (induced by these sections) of 
W we will get such neighborhoods. 

4.1. Bisections of a Lie groupoid. The results contained in this section are well- 
known in the literature (see, for instance, [31 [12]). However, to make the paper more 
self-contained, we will include the proofs of them. 

Let G =^ M be a Lie groupoid with source a : G — t- M and target /3 : G — ?■ M. 

Definition 4.1. A bisection of G is a closed embedded suhmanifold H of G such 
that the restrictions of both a and /3 to are difjeomorphisms . 
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A bisection defines both a section Sq, of a and a section of /3 as follows. 

Proposition 4.2. Given a bisection S the map = is a smooth section 

of a such that ImiJ^a) = S and /3 o is a diffeomorphism. Alternatively, the 
map S/3 = is a smooth section of f3 such that Im(S^) = S and a o is a 

diffeomorphism. 

Proof. Indeed, since S — )■ M is a diffeomorphism, the map is well defined 
and satisfies a o Sq = id-M- Moreover ImE^ = Im((aE)~^) = S and /3 o = 
(/9|s) o is a diffeomorphism because it is a composition of diffeomorphisms. 

The proof of the second statement is similar. □ 
Notice that the diffeomorphisms (3 oTi^ and a o are each one the inverse of the 
other 

{(3 o Sq)""*^ = a o and {a o S/?)""^ = /3 o Sq. 

Definition 4.3. ^4 ZomZ bisection of G is a closed embedded submanifold W of 
G such that there exist open subsets U,V G M for which both a\)^: W ^ U and 
W — V are diffeomorphisms. 

Proposition 4.4. Given a local bisection W the map Wq, = is a smooth local 

section of a defined on the open set U such that (Wq)(W) = W and /3 o Wq : W — )■ V 
is a diffeomorphism. Alternatively, the map = is a smooth local section 

of (3 defined on the open set V such that (>V/3)(V) = W and a o Wp: V ^ U is a 
diffeomorphism. 

Proof. The proof is a straightforward modification of the proof for global bisections. 

□ 

We will need the following straightforward result. 

Lemma 4.5. Let A and B be linear subspaces of a finite dimensional vector space V , 
with dim(A) = dim(i?). There exists a linear subspace G G V such that A (B G = V 
and B®G = V. 

Proof. Let {q} be a basis of AdB. We complete to a basis {q, aa} of A, and we also 
complete to a basis {q, ha} of B. Then {cj, aa, ba} is a basis of A + B, which can be 
completed to a basis {q, Cq,, ba, dj} of V. Then, the subspace G = spanjaa + ba, dj} 
is such that A ® G = V and B ® G = V. D 

Proposition 4.6 (Existence of local bisections). Given g E G there exists a local 
bisection W such that g G W. 

Proof. Since the dimensions of ker(Tga) and ker(Tg/3) are equal, there exists a 
subspace / C TgG such that TgG = ker(Tga)©/ and TgG = ker(Tg/3)©/. Notice that, 
since a and f3 are submersions, we have that Tga{I) = T^^g-^M and TgP{I) = Tf^i^g-^M. 

Let S C G be any submanifold such that (7 G S and T^S = I. The maps T^ajs and 
Tg/3|s are linear isomorphisms at the point g. Indeed Im(Tga|E) = Tga{I) = Ta(g)M, 
and similarly lm(Tgf3\j:) = Tgf3{I) = Tpi^g-^M, and dim(M) = dim(J). By the inverse 
function theorem, it follows that there exist open subsets W" and in S and 
and in M such that g eWnWp and a|w" : >V" -> W° and /3|w/3 : ^ are 
diffeomorphisms. By taking W = H W^, W = a(>V) and V = /3(>V) we have that 
W is a local bisection and g G W. □ 
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In what follows we do not distinguish, in the notation, global and local bisections; 
all then will be denoted by S. 

Definition 4.7. Given a local bisection S defined on the open sets U and V , the local 
left translation by S is the map Ls : a~^(V) — a^^(p() defined by 

Lj]{g) = kg, where h = T,fB{a{g)), 

and the local right translation by S is the map Rj^: fi^^iU) — ?■ fi^^iV) defined by 

R^ig) = gh, where h = Sc,(/3(5')). 

Alternatively, the left action is ■ g = kg, where /i G S is the uniquely defined 
element such that f3{h) = a{g). Similarly the right action is g ■ Tj = gh, where /i G S 
is the uniquely defined element such that a{h) = /3{g). Observe that both Ls and 
i?s are diffeomorphisms. 

It is easy to see that, for a bisection S, the left action preserves the /3- fibers 
and maps a-fibers to ct-fibers. Similarly, the right action i?^ preserves the a-fibers 
and maps /3-fibers to /3-fibers. 

The left action by a bisection S extends the natural left action in the groupoid 
and we have (T, ■ g)h = S ■ (gh), or in other words L-^olg = lT.-g- Moreover, since it 
preserves the a-fibers, it maps left-invariant vector fields to left invariant vector fields 

(35) TLj,(x{g)) = "XiL^g) for every X G r(r). 
Similarly, we have {kg) ■ S = h{g ■ S) or i?^ ° = Rg T.-, from where 

(36) TR,,{l{g)) = t{R^g) for every X G T{r). 

4.2. General discrete Euler-Lagrange evolution operators. Let : G ^ 

be a discrete Lagrangian function and consider a solution (^fQ, Hq) G G2 of the discrete 
Euler-Lagrange equations, that is, 

^(^o)(^d) - ^{h,){L,) = 0, for all X G r(r). 
We denote by xq G M the point 

xq = (3{go) = a{ho). 
Then, we consider the following objects: 

• A symmetric neighborhood W associated with xq and some open subset U of 
G, with local coordinates (x, u) as in Section 13.21 We will denote by V the 
corresponding open subset of M and by (y) the local coordinates on V. 

• A local bisection So of G such that go & Tjq. 

• A local bisection To of G such that ho G To. 

We will assume, without the loss of generality, that the section (So)/3 (respectively, 
(To)a) is defined on the open subset V of M. 

It is clear that Ws,, = -^So(VV) is an open neighborhood of go in G diffeomorphic 
to W. On we may consider a local coordinate system (x, u) defined as follows: if 
g G Wso then there exists a unique gy\; G W such that //^□((y'w) = g] the coordinates 
of the point g are the coordinates in the symmetric neighborhood W of he point gy^. 

Similarly, Wto = RtoO^) is an open neighborhood of ho in G diffeomorphic to 
W. On Wtq we consider local coordinates {y, v) defined as follows: for h G Wtq we 
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consider the unique point hy\> G W such that i?To(^w) = h, and we assign to h the 
coordinates {y, v) of hy^ in the symmetric neighborhood W. 

Moreover, using the same notation as in Section ^72\ the pair of elements {g, h) = 
(Lso((y'w)! -Rto(^w)) £ VVso X Wtq is composable if and only ii y = h{x,u). 

To find the local equations satisfied by the coordinates of a solution [g, h) of the 
discrete Euler-Lagrange equations for La 

(37) t^{g){Ld)-^{h){Ld) = 0, for all 7, 

we take into account equations fl35|) and fl36l) . from where we get that equations f l37|l 
hold if and only if 

iTg^L^o)i^i9w))iLd) - {Th^Rro)ie^ihw))iLd) = 0, for all 7 

or, equivalently, 

^{9w){Ld o ^So) - e^(/iw)(^<i o ^To) = 0, for all 7. 
In conclusion, using f l2^ and f l25|) . we have proved 

Theorem 4.8. The pair {g,h) = (Leo(5'w)) -Rto(^w)) ^ VVe„ x Wt(, «s a solution of 
the discrete Euler-Lagrange equations for if and only if the local coordinates of g 
and h 

g = {x, u) =gw, h = {y, v) = (b(x, u),v) = hw 
satisfy the equations 

(38) LlM?^^^(.,u)+p;(y)?^^i^M+ 



d{LdoRy^) 

-Rliy^^) — ^^-^ — {y.v) = Q. 

Note that the coordinates of go and Hq in and Wtq, respectively, are (a;o,0). 
Next, as in Section [375| we will consider the matrix {¥Ld)J^{x,u), where 

(FL.);:(x,.) = p;(x)^!^^(x,.)+ 



duTi du 



d\LdoRT,) 



Then, we have the following result: 

Theorem 4.9. The following statements are equivalent. 

• The matrix (FLrf)^(xo,0) is regular. 

• The Poincare-Cartan 2-section Q^^ is non-degenerate at the point Hq. 

• The map ¥~Ld is a local diffeomorphism at Hq. 

Any of them implies the following: there exist open neighborhoods Wsq ^ VVsq cind 
Wto ^ VVto of go and ho such that if g & then there is a unique "^{g) = h ^ Wtq 
satisfying that the pair {g, h) is a solution of the Euler-Lagrange equations for Ld. In 
fact, the map : Wsq VVto is a local discrete Euler-Lagrange evolution operator. 
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Proof. If the matrix (FL(i)^(xo, 0) is regular then, using Theorem 14.81 and the im- 
phcit function theorem, we deduce the result about the existence of the local discrete 
Euler-Lagrange evolution operator \E' : — > Wtq- 

On the other hand, the map -Rtq is a diffeomorphism from an open neighborhood 
of e{xo) in the a-fiber a~^{xo) to an open neighborhood of ho in the a-fiber a~^{xo). 
Indeed, (a) it is well defined: ii h E a~^{xo) then 

a{Rro{h)) = a{h ■ To) = a{h) = xq, 

so that i?Xo(^) £ <^~^{xo)', (b) it is injective: it is the restriction of a local diffeo- 
morphism to an a-fiber; and (c) it is surjective: /3 o (To)q, is a diffeomorphism, so 
that if h' G a~^{xo) there exists x E M such that (/3 o (To)a)(a;) = /3{h'), from where 
h = /i'[(To)q(x)]~^ G a~^{xo) and hence -Rto(^) = h'. 
Consequently, if 

d 

w^{ho) = {T,^^o)Rro){j—-^ ^ ), for all 7, 

OU ' \e{xo) 

then {w^lhg)} is a basis of the vertical bundle to a at the point Hq. 
Moreover, if is the right-invariant vector field on Wto given by 

4{h) = {Th^Rr,){4ihw)), for h = Rr,{hw) e Wto 
then, using (1231) . it follows that 

e o _ -p^ — + R^ — 

which implies that the matrix {w^{hQ)(e^{Ld))) is, up to the sign, (FLrf)^(a;o, 0). 

In addition, from ( 1TT|) and (136|) . we deduce that the local expression of the restric- 
tion of ¥~Ld to Wto is 

These facts prove the result. □ 
Now, as in Section [375| we consider local coordinates {x,u) on i{U) = U given by 
(x, u) = (x, u) o i. Note that W C W and, thus, we have local coordinates on Wsp 
and Wto which we also denote by {x, u). 

As above, using equations fl3T]) and fl32|) . we may prove that a pair {g,h) = 
(L^oigw) y Rroi^w)) ^ ^ is a solution of the discrete Euler-Lagrange equa- 
tions for if and only if the local coordinates of g and h 

9 = (2/, v) = (b(S, u),v)= gn>, h = (x, u) = hyv 
satisfy the equations 

(39) = -p'^{bix,u) f^\l^''''\ b{x,u),v) + 



(F"L,)(x,«) = {x,-pl^{x f' '^J^^^\ x,u) + R'^{x,u f' ; ,;^°^ (x,^)). 



dx^ 



+ i?;i(b(x, u),v) ' ^°^ (b(x, ^), ^;) + 
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In order to apply the implicit function theorem (using f l39p . we want to obtain v in 

T 

terms of x and m in a neighborhood of (xq, Uq, Vq)), we consider the matrix (FL)^(a;, m), 
where 

,,.d\LdoRj 



(40) W?S^.u) = ff(x) \LJ"\ ^^^) + 



(x, u) — X, u) + 



Then, one may prove the following result. 



Theorem 4.10. The following statements are equivalent. 



• The matrices (¥Ld)'J^{xo,0) and (FL)^(a;o,0) are regular. 

• The Poincare-Cartan 2-section is non- degenerate at the points Qq and Hq. 

• The maps F+L^ and F~Lrf are local diffeomorphism at and Hq. 

Any of them implies the following: there exist open neighborhoods C VVsq (in-d 
Wto ^ VVx,j of qq and Hq and a unique (local) discrete Euler-Lagrange evolution 
operator \1/ : Wsq VVxo such that '^{go) = h^. In addition, is a diffeomorphism. 



Proof. Suppose that the matrices (FL)^(a;o,0) and (FL)^(xo,0) are regular. Then, 

the existence of the local discrete Euler-Lagrange evolution operator \E' : Wtq 
is guaranteed by Theorem 14.91 Moreover, using fl39l) . fj^Ol) and the implicit function 
theorem, we deduce that there exist open neighborhoods of go and ho which we will 
assume, without the loss of generality, that are Wsq and Wt,, and, in addition, there 
exists a smooth map 

such that, for each h G Wtq, the pair {E{h),h) is a solution of the discrete Euler- 
Lagrange equations for Ld- 

Thus, from Theorem 14. 9[ we obtain that 

(41) oE = id, Eo"^ = id, 

which implies that \1/ is a diffeomorphism. Indeed, ii g & Ws(, then there is a unique 
element h in Wtq, namely h = "^{g), such that {g,h) is a solution of the discrete 
Euler-Lagrange equations for Ld- This proves (HTi) and thus \l/ is a diffeomorphism. 

On the other hand, the map L^.^ is a diffeomorphism from an open neighborhood 
of e{xo) in f3~^{xo) on an open neighborhood of go in /3~^(xo). Using this fact, (fT2|) . 
([21]), (l35l) and proceeding as in the proof of Theorem 14. 9 [ we deduce the result. □ 



4.3. Application. Let L,^ : G — t- M be a discrete regular Lagrangian. Starting 
with a symmetric neighborhood W of the Lie groupoid G, we have local coordinates 
expressions for the discrete Euler-Lagrange equations assuming that we can solve it 
on the symmetric neighborhood. If this is not the case, we are forced to use bisections. 

To illustrate this, assume that G is a connected Lie group and W is a symmetric 
neighborhood of the identity e G G. Assume that there exists a family of points of 
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G: {gi}i<i<m and {^j}i<j<n, with gi = hi = t such that G = [Js'jVV = \^ y^hj (for 

i=i j=i 

instance, if G is compact this property is verified). 

Given an initial point g, then there exists at least an integer I, 1 < I < m such 
that g G giW {gi G G is the bisection using our notation). Now, we try to find 
an integer J, 1 < J < n such that there exits a solution h G Whj of the following 
equation defined on the symmetric neighborhood W: 

iT^-.^L,^){t;igY'g))iL,) - (T,,-.i?,,)(4(/i/i7^))(L,) = 0, for all 7 

where {e^} is a basis of the Lie algebra q. If we find this J, 1 < J < n, we will say 
that the pair {g, h) E G2 is a. solution of the discrete Euler-Lagrange equations for 

We will explore in a future paper these techniques working with points not included 
in symmetric neighborhoods or points not included in the neighborhoods where the 
retraction maps are local diffeomorphisms. 

5. Conclusions and future work 

In this paper we have studied the local description of discrete Mechanics. In Section 
El we have found a local description of the discrete Euler-Lagrange equations using 
the notion of symmetric neighborhood on Lie groupoids. In Section HJ we extend this 
construction for points outside of this type of neighborhoods using bisections. We 
expect that our results apply to a wide range of numerical methods using discrete 
variational calculus. 

On the other hand, this paper will also open the possibility to easily adapt our con- 
struction to other families of geometric integrators derived from discrete Mechanics, 
as for instance, forced or dissipative systems, holonomic constraints, explicitly time- 
dependent systems [U], frictional contact [TB] nonholonomic constraints [7], multi- 
symplectic field theories [13] , discrete optimal control [HI [HI [IZ] . 
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